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Elastic ed scattering, with deuteron polarization, up to high momentum transfer pro- 
vides detailed information on the deuteron wave function. This determines the range 
dependence of the orbital and spin components of the one- and two-body currents, re- 
stricting contributions of isobar and meson-exchange currents and of quark/gluon degrees 
of freedom, as well as the nucleon component. The R- matrix boundary condition model 
combines all these effects, predicting nucleon-nucleon reactions and the deuteron form- 
factors simultaneously. A brief description of the model is followed by a comparison 
of its results with data, emphasizing the restrictions placed on the model by ed elastic 
form-factors. 

1. INTRODUCTION 

There is now elastic electron-deuteron scattering data determining the electric form- 
factor A(q 2 ) up to 6 (GeV/c) 2 U, the magnetic form-factor B(q 2 ) up to 2.8 (GeV/c) 2 [|, 
and the tensor-polarization form-factor t2o(<? 2 ) up to 1.8 (GeV/c) 2 ||. These data restrict 
orbital and spin components of the deuteron wave function at scales as small as 0.2 fm, at 
which distance quark degrees-of-freedom (d.o.f.), isobar components and meson-exchange 
currents all have a significant role. 

The _R-matrix boundary condition method||,[5| provides a hybrid quark/gluon and 
hadron model, incorporating all the above contributions. Only a few of the parameters are 
not predetermined by data independent of the nucleon-nucleon (NN) interaction and sym- 
metry requirements. The remaining few are almost all determined by NN scattering data. 
Essentially one parameter is free to determine the behavior of the three independent elas- 
tic electron-deuteron form-factors (edff) over the large range of momentum-transfers, q. 
This parameter determines the relative amount of AA( 7 Di) and AA( 3 Di) in the deuteron, 
which profoundly affects the q dependence of the spin and convective currents^. The 
NN scattering is not sensitive to the ratio, but only to the sum. 

Following a review of the i?-matrix method and its application to the NN system, three 
specific models for the 1 = 0, J p = I + sector, of different levels of completeness, will be 
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compared with the NN scattering and edff data. From these one can extrapolate to a 
model that represents all the data. 

2. THE i?-MATRIX BOUNDARY CONDITION MODEL 

At high momentum-transfer (short range) the running coupling constant of QCD is 
small, permitting a perturbative description in terms of current quarks and gluons (asymp- 
totic freedom). At low momentum-transfer (long range) nonperturbative effects produce 
clustering into color singlet hadrons (confinement). The transition between these ex- 
tremes has been shown to occur over a small range of the running coupling constant |7[], 
and therefore over a short distance. The .R-matrix method|| is well suited to this situa- 
tion in which two regions, each well represented by a different approximate Hamiltonian, 
have their wave functions connected by a boundary condition at the separating surface. 

For the QCD application, in which confinement requires the quark wave function to be 
small at the transition boundary, the suitable form of the f?-matrix equation is[[|,|| 



Or 

with 



i % 

in which tf) a is the exterior wave function for hadron-pair channel a, W is the total energy, 
the poles Wi are the energies of a complete set of internal states vanishing at r , and the 
residues p l a/3 are given by absolute squares of the internal wave function derivatives at 
the boundaries. Thus the f a p(W) are meromorphic functions with real poles of positive 
residue. The residues can be expressed as 

dWi ■ ■ 

A* = -ro-^Gfi, (3) 

where the fractional parentage coefficients are geometric coefficients expressed in terms 
of Clebsch- Gordon coefficients of the spin/flavor /color space of the given quark configu- 
ration. Therefore only the /°o, representing the effective constant of distant poles and 
the pole at infinity, are free parameters. 

The hadrons at r > r interact via hadron exchange potentials, as given by known 
hadron masses and coupling constants fixed by independent experiments or symmetry 
conditions. 

The separation radius, ro, must be within the range of asymptotic freedom (< 0.85 
of the equilibrium radius of the interior bag model) because of the sensitivity of p l aj3 to 
the derivatives at ro of the internal wave function. Eq. ([!]) also requires that ro satisfy 
^2(^0, Wi(r )) = 0. Using the low energy data, the latter condition fixes the value of r 
with some precision, giving a value consistent with the first condition for the Cloudy bag 
model, but not for the MIT bag model |U|, ruling out the latter for the multi-hadron 
domain. The Cloudy bag model determines r = 1.05 fm. The model then gives a good 
detailed fit to NN data for T La b < 0.8 GeV|J and also is consistent with some evidence for 
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the lowest 1=1 exotic resonance, J p = + |T(| at 2.70 GeV. These resonances, produced 
near the /-poles at Wj correspond to multi-quark configurations other than the minimal 
qq and q 3 configurations. The lowest in the NN system is the I = 0, J p = 1 + at 2.63 GeV. 

3. DEUTERON PREDICTIONS 

3.1. Interior Wave Function 

The first /-pole in the NN system, corresponding to the [g(lsi)] 6 configuration, is in 
the / = 0, J p = 1 + state, 0.76 GeV above the deuteron mass. As the width of the exotic 



resonance is only 0.03 GeV, the f a p are nearly constant. It has been shown |Tl~|| that the 
interior wave function vanishes for constant /, and the actual probability of being in the 
interior has been estimated to be < 0.0040. This implies a large "hole" in the deuteron 
wave function, which has long been noted as a property of the edff and is also embodied 
in the effective "hard core" of NN scattering. In our model this apparent repulsion arises 
from the rapid change in the effective d.o.f. at tq. However at the energy Wi there is 
"matching" of the interior and exterior wave functions, resulting in a substantial interior 
probability and a resonance. 

The fact that r is 2-3 times larger than the cores of repulsive core models is com- 
pensated by the discontinuous increase of wave functions at r , so that the experimental 
effective range parameters can be produced by both types of core effect. However the 
models differ at higher NN scattering energies for large q edff. 

3.2. Exterior Wave Function 

In the I = 0, J p = 1 + system, the NN ( 3 5i) and NN ( 3 L>i) states are coupled to 
each other and to isobar channels by meson exchange potentials as well as the /-matrix. 
Because of their low threshold mass and strong tensor coupling, the AA( 3 Si), AA^D^ 
and AA( 3 D\) channels are most important and are included in all our models. NN 
( 3 <Si) and AA( 3 5i) states have nonvanishing £ l a with the [g(lsi)] 6 quark configuration, 
contributing to the /-pole residue. The NN*(1440)( 3 Si) channel also has a low threshold, 
and with its large width is of next importance to the deuteron. But over the energy 
range which includes the first exotics, other channels are also of some importance and the 
V 1 S n (1535)( 1 Fi), ^ n (1650)( 1 Pi) and AS 31 (1620)( 3 Pi) are included in our recent work. 
These channels modify the best choice of f®» for the lower threshold channels, but have 
negligible components in the deuteron. 

3.3. Determining the /?g for Three Models 

The /9g are sharply restricted by fitting the NN scattering data, which in the deuteron 
sector requires a fit to the np( 3 Si — 3 Di) phase parameters 8 and r]( 3 Si), 8 and ??( 3 -Di) 



and ei for T^b < 0.8 GeV. For the models C and D'|12], which have only NN and AA 
channels, the lower energy behavior of the <fs and e\ determine the NN sector f a p, while 
the energy dependence for 0.4 GeV < TLab < 0.8 GeV fixes those coupling the NN and 
AA( 3 S'i) sectors, and the NN to the sum of AA( 3 A) and AA( 7 L>i) sectors. The last 
two have the same threshold behavior, affecting the elastic scattering in the same way. 
Only detailed A-production data could separate them, so the ratio is free to adjust to the 
edff. The magnetic form factor g-dependence is very sensitive to the ratio because of the 
opposite spin and convection currents of these statesf|. 
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The model C did not consider quark configurations |13| . Without /-poles the best 
choice of separation radius was r = 0.74 fm. Model D' included the lowest /-pole with 
r = 1.05 fm as discussed above. As shown previously 10, C and D' give equivalent fits 
to the 5's while case D' is a better fit to e\. 

When the other isobar channels 
are included, the iViV* (1440), be- 
cause of its larger width, modi- 
fies the energy dependence and in- 
creases the inelasticity. This re- 
duces the required coupling to the 
AA channels. This model E results 
in a better fit to ei,S( 3 Di) and to 
_ 10 the r/'s (Fig. 1). 
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Figure 1. The phase parameters for 
I = 0, J p = 1 + , np scattering. 
Model E (solid curves); SAID SPOO 
phase parameters (dashed curves); 
Bugg 1990-1991 phase parameters 
(squares). 



3.4. The EDFF Predictions 

In previous work[|12|], the edff for models C and D' were calculated with the nonrel- 
ativistic, coupled channel impulse approximation (IA) and the meson-exchange current 
(MEC) terms of it, p, and u, "pair" corrections and the pnj term to first relativistic order. 
For the IA the isobar form factors are assumed proportional to the nucleon electromag- 
netic form-factors. In all cases the MEC corrections to the isobar channels are neglected. 
Both Holder et al. |P.4j| (HO) and Gari-Krumpelmann[|l5f (GK) nucleon form factors were 
used. The results (Figs. 3-6 of flT^fl ) are seen to be a good to A(q 2 ) for the HO choice 
and to B(q 2 ) for the GK choice. This is not inconsistent as A(q 2 ) is dominated by the 
nucleon electric form-factor and B(q 2 ) by the nucleon magnetic form-factor. The t2o(q 2 ) 
predictions were consistent with the very low q experimental results available at the time. 

Here we present, versus the extended data range of the edff, the results of models 
C , D' and E where the first order relativistic correction has been added to the impulse 
approximation and the second order relativistic corrections have been included in the 
MEC[||,0. For model E the ratio of AA( 7 £>i) to AA( 3 £>i) coupling to the NN sector 
was guided by the C and D' model ratios. It has not yet been optimized to the data. 
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Figure 2. A(q 2 ): Data points are described in Ref.[Q. Model C (HO) (solid line); model 
C (GK) (dash-dash); model D' (HO) (dash-dot); model D' (GK) (dot-dot); model _E 
(HO) (long dashes); model E (GK) (dash- dot- dot). 

Figure 3. t2o(o0 : Data points described in Ref. 0. Model C (solid line); model D' 
(dash-dot); model E (long dashes). The dependence on nucleon emff (HO or GK) is 
negligible. 



Also the balance of AA and iViV*(1440) coupling to NN and the value of the iV*(1440) 
magnetic moment, unknown from independent data, have not been varied. 

Table |l] shows the results of model E for the static properties of the deuteron. For 
models C and D' the results are as stated in [12] except for a small relativistic change in 

Q&cut • 



Table 1 

Static deuteron properties of Model E 
BE(MeV) P D (%) P A5 (%) 



Pas(%) Pai(%) Pn*(%T W™ 2 ) Mfr) 



Model E 
Exp t 



2.2247 
2.2246 



5.21 



.006 



3.24 



1.79 



0.71 



.273 
.286 



.860 
.857 



^The higher mass channels have neglibible probability. 



A(q 2 ) is shown in Fig. 2. It is seen that model E with either choice of nucleon form- 
factors fits the data reasonably well for q 2 < 2.5 (GeV/c) 2 , but is only large enough for 
larger q 2 when the GK choice is made. Models C and D' on the other hand are better 
with the HO choice for q 2 < 2.5 (GeV/c 2 ), but at 6 (GeV/c 2 ) only model C with GK is 
large enough. 

^2o(q 2 ) is shown in Fig. 3. For the momentum transfer range there is negligible sensi- 
tivity to the choice of nucleon form-factors, as these cancel in the ratio of quadrupole to 
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monopole electric amplitudes which dominates £20- The result is however very sensitive to 
the model used. The simple constant /-matrix model, C, gives a good fit over the whole 
range of q. Model D' puts the maximum of t 2 Q at much too small a momentum transfer. 
This is related to the large amplitude ofthe L = 2, A A states in this model. Model E, 
with intermediate AA components, has the maximum of £20 between that of models C 
and D'. 

The B(q 2 ) for models C and D' is similar to that of [12] with minima at slightly 



smaller q. However the minimum of B(q 2 ) for model E is at much too small a value 
(q 2 = 1.3 (GeV/c) 2 ). 



4. CONCLUSIONS 

The i?-matrix boundary condition model E, with all the relevant isobar channels, repro- 
duces the np( 3 Si — 3 D\) scattering phases well up to TLab < 0.8 GeV. It also reproduces 
very well the static properties of the deuteron and A(q 2 ) for q 2 < 6 (GeV/c) 2 . It does 
not fit t 2 o(q 2 ) as well as the simpler model C (although it is better than model D'), and 
has the first minimum of B(q 2 ) at much too small a value. To correct the full model (E) 
for B(q 2 ), the ratio of the AA( 7 Di) to AA( 3 Di) couplings to the NN sector needs to be 
varied. That, and perhaps a further substitution of 7V*(1440) coupling for AA coupling 
to the NN sector may also correct the fit to the maximum of t 2 o{q 2 )- 
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